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We model a tiny heat engine as a Brownian particle that moves in a viscous medium in a sawtooth 
potential (with or without load) assisted by alternately placed hot and cold heat baths along its path. 
We find closed form expression for the steady state current as a function of the model parameters. 
This enables us to deal with the energetics of the model and evaluate either its efficiency or its 
coefficient of performance depending upon whether the model functions either as a heat engine or 
as a refrigerator, respectively. We also study the way current changes with changes in parameters of 
interest. When we plot the phase diagrams showing the way the model operates, we not only find 
regions where it functions as a heat engine and as a refrigerator but we also identify a region where 
the model functions as neither of the two. 

PACS numbers: 5.40. Jc Brownian motion-05.60.-k Transport processes -05. 70. -a Thermodynamics 



I. INTRODUCTION 

Understanding the physical properties of devices of 
micron- and nano-meters sizes is of interest these days 
since the trend is in miniaturizing them. Such devices 
may need to be kept under certain environmental condi- 
tion such as a particular temperature in order for them 
to optimally operate. We envision heat pumps (or re- 
frigerators) of similar sizes that provide such favorable 
condition. On the other hand, heat engines of similar 
sizes may provide energy to these devices in order for 
them to perform different tasks. This paper deals with a 
model of such tiny heat engine or, conversely, heat pump 
and explores the details. 

The idea of microscopic (or Brownian) heat engine 
working due to non-uniform temperature first came up 
with the works of Biittiker pj , van kampen 0] , and Lan- 
duaer Q while they were involved in exposing the sig- 
nificance of the now influential papers of Landauer on 
blowtorch effect 0, Q • Later works took models of Brow- 
nian heat engine and dealt with the energetics only at 
the quasi-static limit [E0- Energetics considerations for 
such engine operating in a finite time were first addressed 
by Derenyi et al. and recently by us These two 

works of us also revealed that the quasi-static limits of 
either the heat engine or the heat pump behave exactly 
like that of a Carnot's heat engine or heat pump. 

In our first paper || , we considered a Brownian parti- 
cle moving through a highly viscous medium in a periodic 
sawtooth potential (with or without load) assisted by the 
thermal kick it gets from alternately placed hot and cold 
heat reservoirs along its path. The heat reservoirs were 
placed in such a way that the whole left side of each saw- 
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tooth from its barrier top is coupled to the hot reservoir 
while the whole right side is coupled to the cold reservoir. 

In the present work, instead of coupling the whole left 
side of each sawtooth to the hot reservoir we let it be 
hot within a narrower but varying width located at a 
position that can also vary as a parameter. This will 
introduce two additional parameters to the model and, 
thereby, address a more general problem. By evaluating 
the way energy exchange takes place between the particle 
and the medium, we explore the conditions under which 
the model works either as a heat engine or as a heat 
pump. One new result is the existence of a third region 
in the parameter space where the model works neither 
as a heat engine nor as a heat pump with rich phase 
diagrams. 

The paper is organized as follows: In section II we will 
consider our model in the absence of external load. In 
this case we will show that the model works only as a heat 
engine and, at quasistatic limit, the efficiency goes to that 
of Carnot efficiency. In section III we will consider our 
model in the presence of external load and explore how 
current, efficiency and coefficient of performance of the 
model vary with change in model parameters of interest. 
We will identify the different operating regions (as a heat 
engine, as a refrigerator or neither of the two) of the 
model in phase diagrams. Section IV deals with summary 
and conclusion. 



II. THE MODEL WITH NO EXTERNAL LOAD 

The model consists of a Brownian particle which moves 
in a viscous medium in a periodic sawtooth potential 
where the background temperature is also periodic but 
piecewise constant as shown in Fig. 1. Here U(x) denotes 
the potential of the sawtooth with a barrier height of Uq. 
One of the maxima of the potential is located at x=0 
and the two minima on the left and right sides of this 
maximum potential are located at x = —L\ and x = L2, 
respectively. Two dimensionless parameters, a and 6, 
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FIG. 1: The plot of periodic sawtooth potential without load. 
The parameter a describes the left side of hot locality from the 
nearest barrier top while S describes the width of hot locality 
both in the units of Li. When a — 1 and 5 = the sawtooth 
potential are coupled to the cold reservoir while when a — 1 
and 5=1 the left part of the sawtooth potential are coupled 
with the hot reservoir. 



specify the position and width of each hot locality such 
that a describes the position of the left side of the hot 
locality from the nearest barrier top while S describes 
the width of the hot locality, both in units of L\. Here 
one should note that a and 6 are limited to take values 
between zero and one: < a < 1 and < 5 < 1. On the 
other hand, in our previous work Q both a and S are 
taken to have a fixed value of one. The quantities Th and 
T c in the figure denote the temperature values of the hot 
and cold localities, respectively. A third parameter, r, 
relates the temperature of the hot locality with the cold 
locality such that Th = T c (l + r). The potential, U(x), 
for our model is a series of identical sawtooth potentials. 
One sawtooth potential around x=0 is described by 



1]. 



if -L\ < x < 0; 
if < x < L 2 ; 



(1) 



and this potential profile repeats itself such that U(x + 
L) = U(x) where L = Ly + L 2 . Note that U(x) is de- 
termined by three parameters: Uq, L\ and L 2 . On the 
other hand, the temperature profile is given by 

T(x) = T C + tT c (Q{x - olL x ) - Q(x - (a + 5)Lx), (2) 

for —L\ < x < L 2 and O is the Heavyside function. This 
temperature profile repeats itself with the same period 
as that of the potential: T(x + L) = T(x). 

The presence of hot and cold regions coupled to the 
different parts of the sawtooth potential determines the 
Brownian particle to be driven unidirectionally and at- 
tain a steady state current, J. A general expression for 
the steady state current of the Brownian particle in any 
periodic potential (with or without load) is derived in 
Appendix A. Using the particular potential and temper- 
ature profiles of Eqs. (1) and (2), we have evaluated the 
expressions for F, G\, G 2 and H in Appendix A, Eqs. 



(A10, A11,A12, and A13), and found an exact closed 
form expression for the steady state current J 



-F 



(3) 



The expressions for F, G\, G 2 , and H are given in 
Appendix B. The drift velocity, v, is then given by 



v = J{L X +L 2 ). 



(4) 



As an approximation, we neglect the energy transfer 
via kinetic energy due to the particle recrossing of the 
boundary between the regions 0, @ and find the energet- 
ics of the particle as it exchanges energy with the heat 
reservoirs. In one cycle the Brownian particle is in con- 
tact with the hot region of width 8L\ while the particle 
is in contact with the cold region of width L 2 + Li(l — S). 
When the particle moves through the hot region, it takes 
an energy Qh which enables it not only to climb up the 
potential of magnitude SUq but also acquire energy 5vjL± 
to overcome the viscous drag force within the hot inter- 
val. Hence during one cycle, the amount of heat flow out 
of the hot reservoir, Qh, is given by 



)h = SU Q + JvLyS. 



(5) 



The net heat flow per cycle from the Brownian particle 
to the cold reservoir, Q c , is by a similar argument given 
by 

Qc = 5U -jv(L 2 + L l {l-5)). (6) 
The difference between Qh and Q c is the useful work: 



W = jv(Li + L 2 ). 



(7) 



This is exactly the amount of work required to trans- 
port the Brownian particle moving with a drift velocity 
v through one cycle in the viscous medium. It is impor- 
tant to emphasize at this point that in the absence of 
external load the model works only as heat engine. Fol- 
lowing the definition of generalized efficiency introduced 
in the efficiency of the engine, rj, is given by 



n 



Qh — Qc 



1v(Li 



SiUo + ^vLx) 



(8) 



One can notice that the magnitude of the current ap- 
proaches to zero when Uq goes to zero. This corresponds 
to the quasistatic limit of the engine; i.e. J — » 0. Evalu- 
ating the efficiency at quasistatic limit, we found that 



lim 77 

U o ^0 



(T h - T c ) 



(9) 



which is exactly equal to the efficiency of a Carnot engine. 

In order for us to explore how current and efficiency 
vary with change in one of the newly introduced param- 
eters a and 5, we first introduce three scaled param- 
eters: scaled length £ = L 2 /L\, scaled barrier height 
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FIG. 2: Plot of j versus 8 for fixed t=1, 1=3, u — 4 and 
a = 1. 




FIG. 3: Plot of ?7 versus 8 for fixed r=l, 1=3, u — 4 and 
a = 1. 



u = Uo/ksTc and scaled current j = J/Jq where Jo = 
ksTchLil; Ub being Boltzmann's constant. 

For fixed values of a = 1, u = 4 and r = 1, the current 
j is plotted as a function of <5 as shown in Fig. 2. When 
8 = 0, the whole sawtooth potential is coupled with cold 
reservoir. Hence there is no net current as can be seen in 
the figure. Increasing the width of the hot locality leads 
to an increase in the current. We also plot the engine 
efficiency, r\, as a function of 8 as shown in Fig. 3. The 
figure shows that 77 increases with increase in 8. 

One can fix the width of the hot locality, 8, and in- 
vestigate the influence of shifting the position of the hot 
locality along one side of the sawtooth potential. Figure 
4 is a plot of j versus a for a given values of 8 = 0.2, 
u = 4 and r = 1. The figure shows that the current j 
increases as a increases which implies that positioning 




FIG. 4: Plot of j versus a for fixed r=l, 
8 = 0.2. 



=3, u = 4 and 



FIG. 5: Plot of 77 versus a for fixed r=l, 1=3, u = 4 and 
8 = 0.2. 
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FIG. 6: Plot of periodic sawtooth potential with load where 
periodic hot and cold locality alternate along the potential. 



the hot locality near the potential minimum produces 
the highest possible value of current than putting it any- 
where else. Likewise, plot of r\ versus a in Fig. 5 shows 
similar feature as in Fig. 4 in the sense that one gets 
highest possible efficiency if the hot locality is put near 
the potential minimum than anywhere else. 

In the next section we will consider our model to also 
have external load. 



III. THE MODEL WITH CONSTANT 
EXTERNAL LOAD 

Consider the model in the presence of a constant ex- 
ternal load / as shown in Fig. 6. The potential pro- 
file will then change from U(x) described by Eq. (1) to 
U(x) + fx. Similar to the previous section, the closed 
form expression for steady state current J L in the pres- 
ence of external load is derived and given by 



J L = 



-F 1 



G\G\ + H L F L ' 



(10) 



The expressions for F L , G\ , G\, and H L are given in 
Appendix C. The drift velocity, v L , is then given by 



J i (i 1 +i 2 ). 



(11) 



Before we begin exploring various properties of the 
model, we introduce a scaled quantity associated with 
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FIG. 7: Plots of j versus 8 for fixed r=l, I = 3 ,u 
A = 0.4 and a = 1. 



the load: A = fL\/T c . We will then have a total of six 
parameters that can be controlled independently: u, I, r, 
<5, a and A. Fig. 7 is a plot of (scaled) current j versus 
S for fixed values of all the other parameters. The figure 
shows that the current monotonously increases with 5 be- 
ing negative for small values of S (hot locality width) and 
then becoming positive beyond a certain value of S. This 
implies that in the interval where the current is negative 
external work is being done on the engine while in the in- 
terval where the current is positive the engine does work 
not only against the viscous medium but also against the 
load. In order to find whether the engine works either 
as a heat engine, or as a refrigerator or other wise, we 
need to find energy exchange between the particle and 
the heat reservoirs. When the model behaves as a heat 
engine, the current takes postive value (see Fig. 7) and 
in one cycle the particle takes heat energy from the 
hot reservoir of amount 

Q\ = S(U + fLi + jv L L x ), (12) 

and gives off an amount of heat 

Qc = SU - [((Ml - 6) + L 2 ))(f + 1V L )\ (13) 

to the cold reservoir. The net work, W, done by the heat 
engine in one cycle is given by the difference between 
and Q L C : 



W=( 1 v L + f)(L 1 +L 2 ). 



(14) 



The generalized efficiency of the heat engine will then be 
given by 

( 1 v L + f)(L 1 +L 2 ) 



- YL 

71 ~ ~ 6(U + fL 1 +^yv L L 1 ) 



(15) 



On the other hand, when the model works as a refrig- 
erator heat flows out of the cold reservoir and enters 
the hot reservoir driven by the external work of amount 
W L = f(L\ + L 2 ) done on the refrigerator per cycle. In 
the region where the model works as a refrigrator, the co- 
efficient of performance (COP) of the refrigerator, P re f , 
is given by 

Qc _ SUo - [((Ml -S)+ L 2 ))(f + >yv L )] 
W L /(Lr + L 2 ) 

(16) 
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FIG. 8: Plots of n versus S for fixed r=l, I = 3 ,u = 4, A = 0.4 
and a = 1. 
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FIG. 9: Plots of Pref versus <5 for fixed r=l, 
A = 0.4 and a = 1. 
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The quasi-static limit of the engine corresponds to the 
case where the current approaches to zero either from the 
heat engine side or from the refrigerator side. We have 
found that this limit is satisfied when 



/ 



StU 



((1 + *)(! + t)-Jt)Li" 



(17) 



This is the boundary demarcating the domain of the op- 
eration of the engine as a refrigerator from that as a heat 
engine. Evaluating the expression for both 77 and P re f as 
we approach this boundary, we analytically found that 
they are exactly equal to that of Carnot efficiency and 
Carnot COP, respectively: 



lim r? 



(Th - T c ) 



and 



lim P, 



ref — 



(T h - T c ) 



(18) 



(19) 



We now explore how the rj and P re f behave as a function 
of 5 fixing all other parameters. Within the region where 
the model works as a heat engine, we find the efficiency to 
monotonously decrease from its maximum value (Carnot 
efficiency) as the width, 6, of the hot locality increases 
( see Fig. 8). On the other hand, within the region 
where the model works as a refrigerator the COP linearly 
increases with increase in S until it attains its maximum 
possible value (Carnot refrigerator) as shown in Fig. 9. 

Before we plot the phase diagrams showing regions 
where the model works as a heat engine, as a refriger- 
ator and as neither of the two we think that it is illustra- 
tive to know how Qh and Q c behave with scaled load, A. 
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FIG. 10: (a) Plot q h versus A for fixed r = 4, I = 3, u = 4, 
a = 1 and 5 = 0.6. (b) Plot q c versus A for fixed r = 4, I = 3, 
u = 4, a — 1 and (5 = 0.6. 
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FIG. 11: The phase diagram in the A — r parameter space 
for fixed 8 = .6, £ — 3, u — 4 and a = 1. 



FIG. 12: The phase diagram in the A — S parameter space 
for fixed r=l, I = 3, u = 4 and a = 1. 



M = 4 and a = 1 . One can see in the phase diagram that 
for r = 4, the model works neither as a heat engine nor 
as a refrigerator when 1.21875 < A < 3.156. In fact 
in this region, q c is negative as can be seen in Fig. 10b. 
Figure 12 shows the three regions in A — 5 parameters 
space in which the model operates as a heat engine, as 
a refrigerator and as neither of the two for fixed 1 = 3, 
a = 1, t = land u = 4. 



We introduce two scaled quantities qh = Qh/ksTc and 
Qc = Qc/ksTc corresponding to Qh and Q c , respectively. 
Figs. 10a and 10b show plots of qh versus A and q c versus 
A, respectively, after fixing the other parameters. Note 
that Fig. 10a shows that qh is always positive while Fig- 
ure 10b shows that q c takes a negative value within some 
interval values of A. The quantity qh — q c is always pos- 
itive. In the region that q c is negative, the model works 
neither as a heat engine nor as a refrigerator since under 
this situation external work is supplying energy to both 
the hot and cold reservoirs. We have found that when 



< A < 



Stu 



((1 + 0(1 +t)-*t) 



(20) 



the model works as a heat engine while the model works 
as a refrigerator when 



Stu 



({l + l)(l + r)-5r) 



< A < 



(1-5) + 1 



i(i + 



(21) 

The model works neither as a heat engine nor as a refrig- 
erator if 



A > 



{l-S) + l 



-i(i + 



(22) 



We plot the three operation regions of the model in 
Fig. 11. The figure demonstrates that in the parameter 
space of A — r, the model works as a heat engine, as a 
refrigerator or as neither of the two for fixed 6 = .6, 1 = 3, 



IV. SUMMARY AND CONCLUSION 

We considered a simple model of a Brownian heat en- 
gine and explored how current, efficiency and perfor- 
mance of a refrigerator behave as model parameters of 
interest vary. In the absence of external load, the model 
works as a heat engine while in the presence of external 
load the model works as a heat engine, as a refrigera- 
tor or neither of the two depending on the values in the 
parameter space. At a quasistatic limit the efficiency as 
well as the COP of the engine goes to that of Carnot 
efficiency and Carnot COP. In this paper we reported an 
extended version of the work [9j . The result we obtained 
shows that unlike the previous work, there are regions 
that the model may work neither as a heat engine nor as 
a refrigerator determined by the width of hot locality. 

In this work, we found a closed expression for the 
steady state current. This enabled us to explore the en- 
ergetics of the Brownian heat engine not only at the qua- 
sistatic limits but also while operating at any finite time. 
This is a clear exposition of the power of having analytic 
expression for the concerned physical quantities. 



APPENDIX A: DERIVATION OF THE STEADY 
STATE CURRENT 

Consider motion of a Brownian particle in any periodic 
potential, U(x), and temperature, T(x), profiles where 
both have the same periodic length of L\ + = L as 
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FIG. 13: Periodic potential with non- homogenous tempera- 
ture profile. 



both sides of Eq. (A4) by exp(<f)(x)), we can write it as 
d 



dx 



[e^ x) T{x)P s (x)} = -jJe^ x) . (A5) 



Integrating Eq. (A5) from — L\ to x and, after some 
algebra, one gets 



P°(x) 



e ~4>( x ) 
~T{x) 



T(-Li)P s (-Li) - 7 J J2 L e^dxj 

(A6) 

Taking x=L 2 in Eq. (A6), we get 



shown in Fig. 13. In the high friction limit, the dynamics 
of the Brownian particle is governed by the Smoluchowski 
equation [Til ]: 

where P(x,t) is the probability density of finding the par- 
ticle at position x at time t, U'(x) — ■4-U and 7 is the co- 



( e 0(i 2 ) _ i)ps(L 2 )T{L 2 ) = -7J / " e^dx. (A7) 

Applying normalizing condition, /_£ P s (x)dx = 1, on 
Eq. (A6) gives us the relation 



1 = P S {L 2 )T{L 2 ) 



. T(x) 



■dx 



J J 



L 2 p -<p(x) i-x 



. T(x) 



(A8) 



efficient of friction of the particle. Boltzman's constant, 
ks, and the mass of the particle are taken to be unity. 
The steady state solution, P s {x), of the Smoluchowski Using Eqs , ( A7 ) and ( A8 ) it is s i mp i e to get the steady 
equation implies a constant current, J, given by 



-[U'{x)P s {x) + ±(T(x)P s (x))} = J. (A2) 
7 ax 



Periodic boundary condition implies 
P s {x + L) = P s (x). 



state current 



where 



J = 



dGz+HF' 



(A3) 



Equation A3 is valid for any periodic potential even in 
the presence of additional constant external load. The 
stationary probability distribution can be then given by 



P s (x) 



B 

T{x) 



exp(—(j)(x)) 



(A4) 



F 
Gi 

G 2 

H 



- ML2) 



Li e -4>{x) 



—L\ T[x) 



-dx, 



e^ x '^{x')d'x, 



(A9) 

(A10) 
(All) 

(A12) 



L 2 e -(j>{x) 
Li T(x) 



ilc I e^ x '^{x')d'x. (A13) 



where B is the normalization constant and 4>(x) — 
/-Li U T(xS) wl th — Li < x < L 2 . After multiplying 

I 

APPENDIX B 

In this Appendix, we will give the expressions for F, G\, G 2 , and H which define the value of the steady state 
current, J, for load free case 



F 
G 1 

G 2 



g 5(7o[^--^] _ ^ 

[l_ e ^— ]+ e ^^[l_e r„ ] + 

L, (o-l)U„ UgS (S-a)Ug (°-l)^0 ^0 ^ <.<*-&)Ug Un 

—e T - T h [1 - e T " ] + — e T <= T * T <= [e~ - 1] 

e r 



(Bl) 



(B2) 



(I-q)Uq , UgS (c-S)Ug -u 

e T - + ~ + T - [1 - e~] (B3) 

H = h + t 2 + h + U + h + t 6 + t 7 + t 8 + t 9 + t w (B4) 

*i = ^[(l-a^ + ^Ie 2 ^-!]], (B5) 

t^^e^le^-m-e^], (B6) 

ts = ^[SL 1 + ^[e=^-l]], (B7) 



*7 







-L ' 1 ^. T=i „ T_ T, \„ T=i 



t 4 = -^[1 - e^^]e^^~ — [ e — i^-l], (B8) 



= ^[e^^ - 1][1 - e-]e— , (B9) 



r/ 2 

^0 



Li7 r , „ T r Xi . -^o(°-g) 

*e = ■ 7 r[( a - i ) i i + -^[ c Tc ( B1 °) 



C/o " ' C/o 



rr2 [e^~-l]e T - +T - ^[e--l], (Bll) 



*8 = r : 2 [e T " -l][e--l]e ^ (B12) 
u o 

*9 = rr2 I e ^~ - 1] [e^ - 1], (B13) 

^0 

tio - ^[2*-^[e*&-l]]. (B14) 



APPENDIX C 



In this Appendix we will give the expressions for F L , G\ , G 2 and H L which define the value of the steady state 
current, J L , for nonzero external load case. 

F L = e T - + T- + r h ~c _ i ) (CI) 

Gf = 1-e T c , H e T <= 1-e T '* H e Tc T h 

[1-e -re ] + — 2 e ^ [ e Te _i] (C2) 

„r 7T' f ,_L 1 r (l-°)(t/Q + /I,i) 7T r Ll (1-<*)(E7q+/I.i) , (C/p + /I,i)i (c-aXUp+fL!) lT h L\ (1 - a) (<7 + /L t ) 

-1 + 77T ° ^ ^ +^^+ ^ 1-e— . (C3) 
(Uo - }L 2 ) 

H L = T 1 +T 2 +T 3 + T 4 + T 5 + Te + T 7 + T s + T 9 + T w (C4) 

T2 = (cTtT+T^O 16 ^ - ], (C6) 

T ^raTj [tti + (ra) |e " (C7) 

T 5 = {u ff h f l )2 [e -™>2»~'> _!][!_ e i^^ ]e ^^ ; (C9) 
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T 6 



[(a-S)L 1 + -———[e ^ -1]], 



(Uo + fLi) 

LiL 2 T c 



(U Q + /Li) ' 

(C/ + /J-l)(l-a) 



-(^0 + /^l) | (Up + fLQS (Up + fL^S (U -JL 2 ) 



(Uo + fLMUo-fLtY 



r e T fe - l][e T <= - l]e 



J >> [e J c - 1] 



(Uo-flvWo + fLi) 
{Uo + fLMUo-flvf 1 



[U Q -fL 2 y (U -fL 2 ) 



L 2 T C [e M2)_ 1]L 



lr (Uq-Jl 2 ) 

l][e— ^-1], 



(CIO) 
(Cll) 
(C12) 
(C13) 
(C14) 



r 
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